First principles investigation of ferroelectricity in epitaxially strained Pb2Ti04 
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The structure and polarization of the as-yet hypothetical Ruddlesden-Popper compound Pb2Ti04 
are investigated within density-functional theory. Zone center phonons of the high-symmetry 
K2NiF4-type reference structure, space group M/mmm, were calculated. At the theoretical ground- 
state lattice constants, there is one unstable infrared-active phonon. This phonon freezes in to give 
the 72mm ferroelectric state. As a function of epitaxial strain, two additional ferroelectric phases 
are found, with space groups /4mm and F2mm at compressive and tensile strains, respectively. 

PACS numbers: 77.84.Dy, 77.55. +f, 81.05. Zx 



The occurrence of ferroelectricity in the ABO3 per- 
ovskite structure has been known since the 1950's.^ Re- 
cently, first-principles density functional methods have 
proved invaluable in elucidating the observed behavior of 
known perovskite oxide ferroelectrics, anti-ferroelectrics, 
and quantum paraelectrics^*^ Examples include the 
alkaline-earth titanates,^'^ the alkali metal tantalates 
and niobatesr^ and the lead-based perovskites^iS There 
is also a growing interest in applying these methods to 
the design of new ferroelectrics based on the perovskite 
structure.™'-^ 

Another path to new materials leads beyond the per- 
ovskite structure. The Ruddlesden-Popper (RP) family 
of compounds are closely related to the perovskites:^^ 
They can be viewed as a stacking AO-terminated ABO3 
perovskite [001] slabs of thickness equal to n primitive 
lattice constants. Adjacent slabs are shifted relative to 
one another along [110] by a/2, giving the homologous 
series An+iBnO^n+i- The n=l structure is shown in 
Fig. [2 

Given the structural similarity to the perovskites, it 
seems surprising that there have been no confirmed cases 
of ferroelectricity in the RP family of compoundsii^ Bulk 
phases of RP titanates have been reported only for some 
members of the Sr„+iTi„03„+i and Ca„+iTi„03„+i se- 
riesii^ It should not be surprising that neither the stron- 
tium nor the calcium RP series of compounds appear 
to display ferroelectricity given that the end members 
{n^oo) are SrTiOa and CaTiOs, respectively, neither of 
which themselves are ferroelectric. Still, the fact that 
many RP titanates are not thermodynamically stable 
does not preclude the possibility that a metastable RP 
ferroelectric phase could be produced by an appropriate 
synthetic process. In order to identify such a system, 
it is convenient to use first-principles density functional 
(DFT) methods, for example, to investigate an as-yet 
hypothetical first member of a RP series (where the ef- 
fects of the structural modification would be most severe) 
whose end member is a perovskite ferroelectric. 

In this paper we show that Pb2Ti04, an n=l RP 
compound based on ferroelectric end- member PbTiOs, is 
a promising candidate for a high-polarization ferroelec- 
tric. Regarding the presence of Pb2Ti04 (or higher n 
compounds) in the bulk phase diagram, reports are con- 



fiicting and no structural information is availablei^*i^iii 

However, as discussed above, it may be possible to pro- 
duce a metastable form using modern epitaxial growth 
techniques. Using first-principles DFT calculations, we 
compute the ground-state structure and polarization, 
finding that Pb2Ti04 is a ferroelectric with a polariza- 
tion comparable to PbTiOs. Furthermore, the direction 
of the polarization can be changed by an applied epitaxial 
strain. 

First-principles DFT calculations were performed 
within the local density approximation as implemented 
in the Vienna ab initio Simulations Package (VASP).^^ A 
plane wave basis set and projector-augmented wave po- 
tentials were employed,^^ The electronic wavefunctions 
were expanded in plane waves up to a kinetic energy cut- 
off of 500 eV. Integrals over the Brillouin zone were ap- 
proximated by sums on a 6 x 6 x 6 P-ccntered mesh of 
fc-points. Polarization was calculated using the modern 
theory of polarization as implemented in VASP^ 

We approach the problem of searching for possible 
ferroelectric states by first calculating the properties of 
the RP high-symmetry reference structure, space group: 
lA/mmm, the structure one expects for the paraelectric 
phase. We performed full optimization of the lattice pa- 
rameters (a=3.857A, c=12.70A) and internal coordinates 
(zpb=0.3507, zo„=0.1562). The residual Hellmann- 
Feynman forces were less than 2 me V/A. Next we cal- 
culated the zone-center phonons of this reference system 




FIG. 1: (color online) Crystal structure of Ruddlesden- 
Popper compound Pb2Ti04, Space group I4/mmm. 
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TABLE I: Crystal structure of ferroelectric Pb2Ti04, Space 
Group: /2mm, a=3.985A, fe=3.826A, c=12.70A. 
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by computing the dynamical matrix at g = using the 
direct method where each ion was moved by approxi- 
mately O.OlA. We then froze in the real-space eigendis- 
placements of selected unstable modes and performed full 
relaxations in the space group determined by the symme- 
try breaking mode. Finally, we compute the polarization. 

For Pb2Ti04 in the lA/mmm high-symmetry reference 
structure, there are three infrared-active (i.r.) modes 
that transform according to the irreducible representa- 
tion and four i.r. modes that transform according 
to Eu- The one-dimensional (1-d) ^211 modes involve 
atomic distortions along [001] while in the 2-d £'„ modes, 
atoms move in the plane perpendicular to [001]. Our cal- 
culations reveal that at the ground-state structural pa- 
rameters, lA/mmm Pb2Ti04 has one phonon with an 
imaginary frequency {uj=^&i cm~^), indicative of an in- 
stability. This unstable phonon is infrared-active of E„ 
symmetry type. Therefore, P62 TiO^ in the RP struc- 
ture does indeed display a ferroelectric instability. The 
real-space eigendisplacements of this unstable ferroelec- 
tric mode consists of Pb and Ti atoms moving against a 
non-rigid oxygen octahedra (with larger displacements of 
the apical oxygens in the PbO layer). As for PbTiOs,^ 
the character of the ferroelectric instability in Pb2Ti04 
involves significant Pb displacements moving against oxy- 
gen in the Pb-0 planes. This involvement of the A-site 
cation in both PbTiOa and Pb2Ti04 differs from non-Pb 
based compounds (e.g. BaTiOs) and has been attributed 
to the Pb^+ 6s^ lone-pair. It is in fact this lone-pair 
physics that stabilizes PbTiOa in the tetragonal phased 
and may have a role in facilitating the ferroelectric dis- 
tortion in Pb2Ti04.2i 

The key role of Pb is further emphasized by comparison 
with Ba2Ti04. A compound at this composition, barium 
orthotitanate, has been identified in the bulk phase dia- 
gram. It crystallizes not in the RP structure but rather in 
the monoclinic distorted-K2S04 structure (space group: 
P2i/n).^^ A calculation for the structure and zone-center 
phonons of lA/mmm Ba2Ti04 in the hypothetical RP 
structure, exactly analogous to that for Pb2Ti04, shows 
no evidence for any ferroelectric instability, even for vary- 
ing epitaxial strain. 

Returning to Pb2Ti04, we search for the ferroelectric 
ground state by freezing-in the real-space eigendisplace- 
ment pattern of the unstable E„ mode, and perform- 
ing full relaxations of all ions and lattice constants con- 



sistent with the resultant space group. Since this E„ 
mode is doubly degenerate, any linear combination of 
the degenerate modes polarized along [100] and [010] is 
an equivalent eigendisplacement. We considered two lin- 
ear combinations; one polarized along [100], a second po- 
larized along [110]. Freezing-in the E^ mode polarized 
along [100] results in a body-centered orthorhombic space 
group, I2mm. For the distortion along [110], the result- 
ing space group is face-centered orthorhombic F2mm. 
Since this F2mm structure is slightly higher in energy 
(20 meV/formula unit) than the I2mm structure, we now 
consider only the latter (we will revist F2mm below when 
discussing the effect of epitaxial strain). Our calculated 
structural parameters of Pb2Ti04 in this orthorhombic 
space group are displayed in Table |l| We imposed the 
convention, Axi=0, where Axi is the deviation along 
[100] from the centrosymmetric position of ion i. It can 
be seen that the relaxed structure can in fact be related 
to the high-symmetry RP structure by small displace- 
ments of Pb and Ti ions moving against the non-rigid 
oxygen octahedra, consistent with the freezing-in of the 
E„ phonon instability of the high-symmetry structure as 
proposed. Finally, we calculate the spontaneous polar- 
ization Ps. We find that P^ = GS^C/cm^, along [100] 
as required by symmetry. Therefore P62 TiOi in the RP 
structure is a ferroelectric with a spontaneous polariza- 
tion comparable to that of PbTiO^?^ 

Epitaxy plays a dual role in our thinking about the 
Pb2Ti04 system. As will be discussed shortly, one possi- 
ble route to synthesize thin films of Pb2Ti04 in the RP 
structure is through the use of epitaxial stabilizationi^iiS^ 
In addition, it is becoming increasingly possible to grow 
oxide thin films coherently on substrates with a relatively 
wide range of lattice constants (1-2% lattice mismatch is 
currently the norm) . This provides an additional param- 
eter to "tune" the properties of the material to desired 
values by applying an in-plane (or epitaxial) strain to the 
thin film compared to bulk. 

With this in mind we consider again the high- 
symmetry, lA/mmm RP structure and explore the effect 
of epitaxial strain on the low-frequency infrared-active 
modes. We impose epitaxial strain by constraining the 
two basal primitive vectors of the bet lattice to an angle 
of 90 degrees and to a fixed equal length (i.e. correspond- 
ing to that of an implicit coherently matched square- 
lattice substrate) . In Fig. |21 we show how the phonon 
frequencies for the lowest frequency E„ and A2M phonons 
vary as a function of compressive epitaxial strain. We 
use the computed a parameter of tetragonal ferroelectric 
PbTiOj^'^ as the reference strain (i.e. for 0% strain we 
fixed the in-plane lattice constant of Pb2Ti04 to that of 
ferroelectric VhT\Oz)t^ For each value of fixed strain we 
again perform relaxation of the ions and c-axis. 

Referred to tetragonal PbTiOs, the ground-state 
lA/mmm structure has an in-plane strain of approxi- 
mately -0.3%. For small compressive strains the phonon 
instabilities of epitaxial Pb2Ti04 are expected to be sim- 
ilar to those in the unconstrained structure. This corre- 
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FIG. 2: Soft infrared-active phonon frequencies as a func- 
tion of in-plane compressive strain for the lowest frequency 
Eu phonon and lowest frequency A2u phonon of the high- 
symmetry, paraelectric lA/mmm reference structure. 



sponds to the far right of the Fig. [21 It is evident that the 
lattice dynamics in this region of zero or small epitaxial 
strains are dominated by the largely unstable E„ mode, 
as previously discussed. If we now increase the compres- 
sive epitaxial strain (i.e. from right-to-left) the E„ mode 
stiffens while the mode softens considerably and be- 
comes unstable at « -2.5% strain. Fig.[21shows that for 
large compressive strains (-4 to -5%), the highly unsta- 
ble should dominate the lattice dynamics while for 
intermediate strain values, both an mode and an E„ 
mode are unstable and comparable in value. This behav- 
ior with strain can be simply understood as arising from 
volume effects. As we increase the in-plane compressive 
strain, the effective volume in which the mode (po- 
larized in-plane, perpendicular to the c-axis) vibrates de- 
creases. This increases the short-range repulsive forces, 
thereby stiffening the force constant i2i In contrast, the 
effective volume of the A2„ mode (polarized parallel to 
the c-axis) increases with increasing compressive strain 
leading to a softening of the force constant. 

Next we use these phonon instabilities (Fig. as a 
guide to search for additional epitaxially stabilized ferro- 
electric structures. At each value of strain we first freeze- 
in separately the real space eigendisplacement pattern 
corresponding to the A2„ mode {I4mm), the E„ [100] 
mode (/2mm), the E„ [110] mode (F2nnn), and both 
the A2„ and the E„ [100] modes {Cm). Then we relax all 
ions and the c-axis lattice parameter while keeping the 
in-plane lattice parameters fixed. Finally we calculate 
the total energy. Fig. ^ and the polarization. Fig. 13 of 
the resultant structures as a function of epitaxial strain. 

As shown in Fig. O RP Pb2Ti04 undergoes a series 
of structural transitions with epitaxial strain. Over the 
range of slightly tensile to compressive strains, Pb2Ti04 
forms in the /2m,m space group. Following the conven- 
tion appearing in the literature^ we refer to this phase 
as the a-phase. The polarization in this phase varies 



FIG. 3: Polarization along [001] {Mmm), [100] (/2mm), and 
[110] {F2mm) as a function of epitaxial strain. 



from 34^C/cm2 at -3.3% strain to SG^iC/cm^ at -1-0.7% 
strain. The minimum energy structure in the a-phase oc- 
curs at a tensile strain of -|-0.55%, corresponding to the 
ground state structure discussed above. As the compres- 
sive strain increases, the energy of the /2mm structure 
approaches that of the paraelectric lA/mmm structure, 
as shown in Fig. ^ while remaining about 3 meV lower 
for the values of strain that we considered. This is con- 
sistent with the leveling off of the E„ phonon as shown 
in Fig.|2]and explains why the polarization of the /2mm 
structure remains nonzero for large compressive strains. 
A transition from the a-phase, with the polarization in- 
plane, to a phase with the polarization along the c-axis, 
i.e. from /2mm to /4mm, occurs for large compressive 
stains as anticipated from the phonon instabilities of the 
I A/ mmm structure. This occurs at Ri-3.3% strain. We 
refer to this lAmm phase as the c-phase. The polarization 
in the c-phase approaches 70/iC/cm^ at -4.0% strain and 
continues to increase. In the range where the two unsta- 
ble mode frequencies cross, we considered the possibility 
that coupling between the two modes could lead to ad- 
ditional ferroelectric structures. However, relaxing the 
structures in the low symmetry space group Cm always 
yielded one of the two higher symmetry structures (a or c 
phase, depending on the value of the misfit strain). Thus 
the transition from a to c with increasingly compressive 
in-plane strain appears to be first-order.^^ Finally, for 
large enough tensile strains (greater than wO.7%) the 
F2mm structure becomes lower in energy than that of 
the /2mm structure. The polarization of this aa-phase 
is comparable to that of the a-phase while the minimum 
energy structure occurs at a slightly more positive strain 
of -1-0.7% strain. The point at which the energy curves for 
the a and aa phases cross is of particular interest, as the 
in-plane polarization is nearly isotropic. The free rota- 
tion of the polarization might result in some unexpected 
interesting physical properties. 

Compounds unstable in the RP structure at atmo- 
spheric pressures have been synthesized under high- 
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FIG. 4: Energy (per formula unit) as a function of epitaxial 
strain for space groups lA/mmm (paraelectric), 72mm (a- 
phase), F2mm (aa-phase), and 1 4mm (c-phase). 

pressure conditions, e.g. polycrystalline Ba2Ru04.'^'' 
This route seems promising to synthesize bulk RP 
Pb2Ti04. Further, we suggest synthesis of non-bulk 
phases of this material through the use of epitaxial sta- 



bilizationiMiS^ In fact, this method has proven quite suc- 
cessful for stabilizing a variety of oxide thin films in the 
RP structure. One example is the higher order mem- 
bers {n > 3) of the Sr„_|_iTi„03,i+i series where bulk 
phases are known to be unstable. '^^ Another example has 
been the low pressure synthesis of Ba2Ru04 in the RP 
structure. •^■^ The growth of thin films of Pb2Ti04 would 
provide a means to realize the interesting behavior of this 
material with epitaxial strain. 

Initially, we asked the question whether Pb2Ti04 in 
the Ruddlesden-Popper structure would be a ferroelec- 
tric. Using first-principles DFT calculations, we have 
seen that indeed it does display a ferroelectric instabil- 
ity. We have argued that if Pb2Ti04 could be made in 
the RP structure (bulk or thin films) it would undergo 
a ferroelectric structural transition to the orthorhombic 
a-phase with a spontaneous polarization comparable to 
that of bulk PbTiOs. 
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